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Abstract 

We measure the torque acting on a rod rotated perpendicular to its axis in a granular bed, 
through which an upflow of gas is utilized to tune the hydrostatic loading between grains. At low 
rotation rates the torque is independent of speed, but scales quadratically with rod-length and 
linearly with depth; the proportionality approaches zero linearly as the upflow of gas is increased 
towards a critical value above which the grains are fluidized. At high rotation rates the torque 
exhibits quadratic rate-dependence and scales as the rod's length to the 4th power. The torque 
has no dependence on either depth or airflow at these higher rates. A model used to describe the 
stopping force experienced by a projectile impacting a granular bed can be shown to predict these 
behaviors for our system's geometry, indicating that the same mechanics dictate both steady-state 
and transient drag forces in granular systems, regardless of geometry or material properties of the 
grains. 
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I. INTRODUCTION 



Granular systems are ubiquitous in both nature and industry. While such systems are 
comprised of many mechanically simple particles, their bulk properties tend to be complex 
and nonlinear 

HQS 

, |4j. For example, a static pile of sand might behave as a solid, but 
if sheared or shaken strongly enough the pile will spontaneously begin to flow much like a 
liquid |5|] . This capacity to flow makes granular systems both interesting and useful, however 
traditional rheo logical models are not sufficient to fully describe granular flow. As a result, 
developing continuum approximations for granular flow is a long-standing goal. 

An intuitive quantity to study in the context of granular flow is the drag exerted by 
a granular bed on a moving intruder. A general characterization of such resistive forces 
could be valuable in numerous applications - such as predicting the outcome of a projectile 
impact into granular objects such as asteroids or planetoids jf], Q, [s,^, 10, 11] or determining 



the forces experienced by industrial mixing and digging equipment - and may contribute 
to a more comprehensive understanding of granular matter [12], Q, Q, H, 1 1 €31 ] . Previous 
attempts to formulate such a general characterization in the context of impact dynamics 
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have resulted in several contradictory models for the drag force jll . 

In a recent study of impact dynamics for a 2.54 cm diameter steel sphere, Katsuragi 
and Durian [18j demonstrated that the stopping force scales as the sum of two terms, one 
which grows linearly with depth, z, and one that grows quadratically with velocity, v. This 
functional form succesfully describes not only the results obtained in Ref. [3], but also the 
observations reported in the numerous impact experiments which had been held to support 
one or another of the aforementioned contradictory drag force models. 

The scaling of the two force terms observed in Ref. [l8| suggests the following prediction 
for the materials dependence of drag force in the bed. The force component proportional 
to v 2 resembles the inertial drag experienced by an object moving through a fluid at high 
Reynolds number. Therefore, this force term ought to scale with the cross-sectional area 
of the intruder, A, and the density of the fluid - in this case the mass density, p g , of the 
granular packing. The force component proportional to z is reminiscent of Coulomb friction. 
Accordingly this force term should scale like a friction coefficient, fi = tan (6 r ) where 9 r is 
the angle of repose, times a hydrostatic pressure, p g gz, where g is acceleration due to gravity. 
Similar results have been obtained for the drag acting on intruders steadily driven through 
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grains in the slow, quasistatic limit 19j, |20_|, |2l|, |22j. By this reasoning we expect the full 
material dependence of the drag force to be of the form 



F = afxp g gAz + Pp g Av 2 



(1) 



In Ref . 



181 ] the numerical prefactors were determined by fit to be a = 26±3 and /3 = 1.0±0.1. 



If the materials-dependence of Eq. flTJ is correct, these values should be independent of the 
details of the system, however this has yet to be tested. 

While impact is a natural context in which to observe granular drag on an intruder, 
there are some drawbacks to using impact experiments to carefully characterize drag forces. 
Impacts are transient, with conditions that are inherently time dependent. Variables such as 
velocity and depth are entangled in an impact, adding uncertainty to the observed scaling of 
stopping force with these parameters. Furthermore, the drag force is generally determined 
by measuring instantaneous position or velocity rather than by direct measurements of 
stopping force, thus limiting the force resolution of experimental data. In order to address 
these limitations we develop a 'steady-state impact' experiment by rotating a horizontal rod 
that is submerged in a granular bed, and measuring the torque acting on the rod once the 
drag has reached a steady-state. This design enables us to vary control parameters such as 
rate and depth in isolation, and to make measurements of drag torque directly. 

Our apparatus is similar to geometries that have been utilized for the study of granular 
flow in the past. The most prevalent such systems are bladed mixers and vane rheometers. 
Bladed mixers consist of short, angled vanes rotated or pushed through a granular bed. 
Studies utilizing bladed mixers tend to focus on the geometry of the flows generated by 
the motion of the blades, emphasizing mixing efficiency rather than the explicit form of 



the torque acting on the blades 
number of vertical vanes, general 
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261 ] . Vane rheometers are devices in which a 



about a vertical axis 
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y spanning the entire depth of a granular bed, are rotated 



301 ] . Vane rheometers attenuate the complications caused 



by wall-slip in the more rheologically conventional Couette geometry, and a re p rimarily 
utilized to characterize yield stresses, thixotropic properties, and flow-profiles 27|. In this 
geometry, rates and torques agree qualitatively with the model presented in Ref. [18], but 
are considered in terms of stresses along the surface of the c ylin der of granular media that 



is imagined to rotate as a solid along with the vanes 
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30J. In contrast, in this paper 



we picture the forces created by the motion of an intruder moving through a medium that 
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is otherwise at rest. 

In addition, we seek to test our interpretation of the mechanical origins of granular drag. 
As discussed above, contact forces in the granular bed are assumed to scale as the hydrostatic 
pressure independent of force chains activated by intruder motion. We modify this pressure 
by generating a flow of gas through the granular bed which, by Darcy's law, generates a 
pressure drop proportional to the gas speed through the bulk. We extend the model to 
account for this effect by making the substitution: 

/ (1 - U/U c ) Pg gz U<U C 

Pg9 z -> < (2) 
unknown U > U c 

where U is the superficial gas speed equal to the gas flux divided by the sample cross- 
section, and U c is a critical gas speed which is sufficient to fluidize the bed. According to 
this substitution the only role the airflow plays is to reduce contact forces, and will thus 
have no effect on the inertial component of the drag. Air has been shown to play a subtle 
but dramatic role in the formation of convection-like patterns in shaken systems [3JJ and the 
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331 ] . so such a simple dependence 



dynamics of granular jets generated after an impact 
on gas speed is an interesting result in and of itself. For U > U c we expect a qualitative 
change in dynamics, and that the the system will cease to have a yield stress. 

We verify that our experimental observations agree with this drag force model. These 
results reinforce the supposition that the dynamics of drag are a universal characteristic of 
granular systems, and that dynamics at low rates are set by static contact forces while for 
rapid perturbations drag is predominantly inertial, and that these dynamical regimes are 
universal traits of granular drag. 



II. MATERIALS AND METHODS 



The experiments were performed with an Anton-Paar UDS200 rheometer. The rheometer 
drives the rotation of a horizontal rod, and measures the torque and corresponding rotation 
rate experienced by the rod. The rheometer provides us with resolution in torque as fine 
as 10~ 5 Nm, and rate resolution of 10~ 4 rotations per second. We are able to access a 
dynamical range of five orders of magnitude in both rate and torque. 

As depicted in Fig. [TJ the rheometer tool is replaced with a horizontal rod which rotates 
about its center of mass, perpendicular to its axis. The data we present were collected using 
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three different rods with diameters (lengths) of 12.6 mm(76 mm), 6.34 mm(102 mm), and 
6.34 mm(96 mm). 
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FIG. 1: This is a sketch of the dish used to apply a pressure drop across the granular bed. The lower 
chamber, or windbox, is pressurized with air, which results in a homogenous flow of air through 
the porous glass frit. During an experiment the upper chamber is filled with spherical glass beads. 
The upflow of air through the grains generates a pressure drop which reduces hydrostatic loading. 
A horizontal rod of length L c and diameter D c is plunged horizontally into the glass beads to depth 
z, then rotated perpendicular to its axis. 

In place of the rheometer's standard sample plate we mount a 12.7 cm diameter dish, 
also depicted in Fig. [TJ The upper portion of the dish is filled with spherical glass beads - 
Potters Industries with radii of 125 fim (stock number P-0060) and 180 [im (P-0080) - to a 
depth of 10 cm. The lower portion is a windbox, 13 cm tall. Separating the two chambers 
is a porous glass frit. The frit is a Robu and Schmidt P2 porous glass disc with a thickness 
of 8mm, and pore-sizes ranging from 50 to 100 (im. A pressure drop across the glass frit 
produces a homogenous upflow of gas into the granular bed. This flow of gas through the 



sample has been shown to produce a uniform pressure gradient which opposes gravity [34j. 
and provides a means to tune the magnitude of static contact forces in the bulk. We use 
air to pressurize the windbox, and the air is passed through a pressure regulator, desiccant 
air-dryer, and adjustable airflow meter before entering the windbox. The pressures in the 



windbox are sufficient to generate superficial gas speeds of 0-44 mm/s, well in excess of 
the critical fluidization airspeeds. The critical airspeed for the onset of fluidization was 
determined to be U c = 13.7 mm/s for 125 /zm grains and U c = 22.5 mm/s for 180 //m grains. 

To prepare a granular bed we begin by fluidizing the grains with dried air before beginning 
the experiment in order to mitigate any ambient humidity. We then set the gas speed through 
the sample and the depth of the rod to values for which we wish to conduct a measurement. 
This apparatus enables us to make two different types of measurement: static measurements 
in which we slowly increase the torque from zero until reaching the yield torque r y at which 
the rod begins to rotate, and dynamic measurements in which we set the rod to rotate 
at a series of rates and record the mean torque required to maintain that rate. When we 
run our experiment at a constant torque, measurements of the instantaneous rotation rate 
exhibit fluctuations of order 10%, which are largest at low rotation rates, and vanish at high 
rotation rates. Since granular systems don't exhibit conventional rheological flows, nor does 
our apparatus have a conventional rheometric geometry, we present our results in terms of 
torques and rotation rates rather than stresses and strain rates. 

III. DATA 

Example data for torque vs rotation rate are presented in Fig. [2} with upper and lower 
plots being for different grain and rod sizes, and where each curve corresponds to different 
airspeeds as labeled. At low airspeeds, and for low rotation rates, the torque is constant, 
hence the behavior is quasistatic. The values of the torque in this limit, r , decrease sys- 
tematically as the upflow of air is increased. They correspond well to the yield torques, r y , 
indicated by the solid triangles at zero rotation rate. At airspeeds above the critical airspeed 
U c , at which the granular bed begins bubbling by eye, r y vanishes, though a nonzero torque 
is measured during rotation that decreases with rotation rate. At very high rotation rates, 
generally above one rotation per second, the measured torque increases rapidly with rotation 
for all airspeeds. 

The behavior in the quasistatic regime is compared with expectation in Fig. [31 where 
both To and r y are plotted against the normalized airspeed U/U c for the same conditions as 
in Fig. [2j As noted already, the values of tq and r y are in fair agreement, decreasing with 
U, and vanishing for U > U c . Furthermore, the dashed lines in Fig. [3] indicate that the 
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FIG. 2: Torque plotted against rotation rate for gas speeds ranging from to 44 mm/s. The 
top [bottom] plot is data from experiments conducted with 180 [125] /xm grains, for which the 
fluidization (boiling) gas speed is around 22.5[13.7] mm/s at a depth of 50 mm. The dashed red 
lines are fits to Eq. ([3]), and are all constrained to have the same value of (3, corresponding to the 
solid black lines. The zero-rate valued triangular points denote the yield torques measured under 
the same conditions as the dynamics data for which fits are shown. 

quasistatic torque approaches zero at fluidization in proportion to (1 — U/U c ). This scaling 
is precisely that expected by Eq. in which grain-grain contacts are loaded hydrostatically 
by gravity with an offset proportional to airspeed. 

To more fully test the dependence of the quasistatic torque on airspeed as well as on 
the details of the system, and to investigate the behavior at higher rotation rates, we have 
collected over one hundred sets of torque data vs rotation rate using 3 rods, 2 grainsizes, 
and 40 depths ranging from 2-8 cm. In all cases the results are qualitatively similar to those 
in Figs. EH 
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FIG. 3: ro and r y plotted against gas speed divided by the critical fluidization gas speed for both 
180 fim (top) and 125 /jm (bottom) diameter grains. The lines are linear fits to the low gas speed 
data. D c is the rod diameter, and Lc is its length. 

IV. ANALYSIS 

To test the force model of Eqs. (|If2]) . and ultimately deduce the values of the proportion- 
ality constants a and (3, we compute the torque as J A rdF and compare with data. In the 
limit that the rotating rod is long and thin, the speed and area of a differential element at 
distance r from the rotation axis are v = lot and dA = D c dr, respectively, where u is the 
angular rotation speed and D c is the rod diameter. Substituting into Eq. (p]) and carrying 
out the integration along the length of the rod gives the following prediction for the total 
torque produced by granular drag: 

r = r + -^p g D c L 4 c u 2 , (3) 

a ( U \ 2 
to = 4 ( 1 ~ ) W 9 gD c L c z, (4) 

where L c is the rod length. Fits to this form are shown as dashed curves in Fig. [2j where 
t was adjusted independently for each data set at U < U c but where a single value of (3 
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was enforced by simultaneous fit to all data for a given grain size and probe geometry. The 
inertial contribution to the torque is plotted as a solid black line. The high quality of the 
fits with a single value of (3 clearly demonstrate that the force model of Eq. (pQ) accurately 
approximates the behavior of the medium. 




(1/4)(1 -U/U> Pg gD c L c 2 z (Nm) 

FIG. 4: ro plotted against (1/4)(1 - U /U c )np g gD c L 2 c z for both 180 ^m(top) and 125 //m(bottom) 
grains. Coloring denotes unsealed gas speed, and symbol size denotes depth, where the largest 
symbols correspond to a depth of 78 mm, and the smallest correspond to depths of 20 mm. 

Next we focus on experimental data for the quasistatic torque, r , and their comparison 
with the explicit materials- and geometry-dependence of Eq. (J3J). For each data set, the 
value of ro is found by fitting the low rotation rate data to a constant. The value is taken 
as ro = for airspeeds where the torque systematically decreases as to is reduced, as for 
U > U c . Results for all runs are plotted in Figs. HK,b, for the two different grain sizes, as a 
function of x — (1/4)(1 — U/U c )iJ,pggD c L c 2 z. According to the prediction of Eq. (@J, both 
axes have units of torque and the results for x < should be zero. But more importantly, 
Eq. (jlj) predicts that all data for x > should collapse to a line through the origin with 
slope a. Indeed this is consistent with observation, showing that Eq. (jlj) correctly predicts 
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the dependence of the quasistatic torque on rod diameter, rod length, depth of rotation, 
and the speed of upfiowing air. Furthermore, the fits to tq = ax shown as lines in both 
Figs. @Ji,b give the same value a = 15 ± 2 for the two different grain types, suggesting that 
Eq. (J3J) captures the complete material dependence of the quasistatic torque as well. 

While the analysis of quasistatic torque thus appears in good agreement with expectation 
for the materials and geometry dependence, two words of caution are in order. First, the 
deduced value a = 15 ± 2 is somewhat smaller than the result a = 26 ± 3 found in Ref. [3] 
based on the dynamics of impact of a single steel sphere into a single medium for a range 
of drop heights. This could reflect actual differences due to the direction and/or geometry 
of intruder, horizontal rods here vs downward spheres in Ref. jl8|; or it could reflect an 
uncertainty in Ref. due to the entangling of the position- and speed-dependent force 
terms inherent in impact. Second, the linear behavior of To vs x does not appear to hold 
near x = 0, and nonzero r values are even found for x slightly less than zero. We believe 
both effects are due to small manufacturing defects in the porous glass air diffuser, such 
that the speed of air is not perfectly uniform across the bottom of the sample. A thinner 
region would give a higher local airspeed, and fluidize a portion of the sample at U < U c ; 
the same effect would occur for a site that persistently nucleated bubbles. Unfortunately 
this is difficult to explore because the granular medium is opaque. 

Lastly we consider the high rotation speed torque data and its comparison with the 
prediction of Eq. ([3]). For this, we simply subtract the quasistatic torque and plot the 
result r — To as a function of x = (l/32)p g D c L c u> 2 in Figs. 5a,b for the two different grain 
types. According to Eq. ([3]), both axes have units of torque and all data should collapse 
to the line r — r = f3x. Indeed this is consistent with observation, showing that Eq. (J3]) 
correctly predicts the materials and geometry dependence of torque in the inertial regime. 
Furthermore, the fits to a line are consistent with a single value j3 = 0.6 ±0.2 for both grain 
types. This is somewhat smaller than the result j3 = 1.0 ± 0.1 of Ref. 18| based on impact, 
perhaps for one of the reasons noted above. Note, however, that since a and (3 are both 
smaller for rotating rods than for downward spheres, it seems fair to conclude that the drag 
on the latter is larger. 



10 




0.05 0.1 0.15 
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FIG. 5: Torque plotted against (1 / '32) pD c L^.uj 2 , and collapsed by subtracting off the quasistatic 
term from Eq. ([3]). The data on the top [bottom] plot is from 180[125] /jm grains, depths between 
20 and 78mm, and gas speeds between and 44mm/ s. For U > U c we take To to be 0. The lines 
have slopes equal to (3 and go through the origin. 

V. CONCLUSION 

To summarize, we have developed a new method to study granular drag. Our system 
achieves a reproducible steady-state, facilitating statistically robust data collection, and 
enables us to vary variables such as speed and depth independently and in situ. This system 
has the potential to be a powerful tool in understanding the fundamental bulk mechanics of 
granular flow: already our results have shed light on the grain-scale mechanics which dictate 
drag forces on the scale of an intruder. We observe drag forces proportional to the rate 
squared for high rotation rates with magnitudes which depend only upon the geometry of 
the intruder and the density of the granular packing, suggesting drag forces at high speeds 
are inertial in origin. In the slow quasistatic limit we confirm that the drag force is rate- 
independent, and scales like a hydrostatic pressure. This scaling suggests that the response 
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of a granular packing to low-speed intruders is set by the strength of static contact forces. We 
modify the loading of contacts between grains in the packing with an applied pressure drop 
across the bed, and the agreement between data and Eq. <^ provides further corroboration 
for our hypothesis. A simple model based on this interpretation of the mechanistic origins of 
granular drag has proven effective in describing impact dynamics [18|. That we have similar 
success in applying the model to describe the drag experienced by an intruder moving 
horizontally at a constant rate reinforces the supposition that the nature of granular drag is 
universal, independent of the details of the system. 

An obvious open question is why the values we obtain for a and 3 are different from those 



reported in [18]. This disparity suggests that Eq. fl3]) might not capture the full geometry 
dependence of the drag force. As discussed above, one such geometric effect might be the 
dependence of drag forces on the direction of intruder motion relative to the direction of 
gravity. Another aspect of geometry that might warrant closer examination is intruder 
shape. Cross-sectional area may not be sufficient to fully describe the role of intruder shape, 
and corrections for three-dimensional features of intruders may prove important. A further 
characterization of granular drag forces in a more extensive set of geometries is needed to 
fully determine the role of such geometric effects. 

Another area of research that remains largely unexplored is the full form of the drag force 
in the fluidized case. While gas-ffuidization is a fairly specific phenomena, it's also a useful 
way to induce flow, and understanding the non-inertial components of the drag force in a 
gas fluidized bed may provide insight into the behaviors of other granular systems in which 
there are no lasting static contacts and grain-grain interactions tend to be ballistic, such as 
in very rapid and low-density granular flows. As such, an in depth investigation of low-rate 
drag forces for U > U c would be of great interest. 

Lastly, understanding the effect of more complicated grain interactions upon the drag 
dynamics in a bed would be essential to describe drag in many real-world granular systems. 
Our model assumes that the only mechanisms that contribute to drag are inertial interactions 
between the intruder and the grains, and static contact forces between grains. For our system 
- mechanically simple, hard, dry, spherical grains - these are good assumptions, but many 
granular systems are not nearly so simple. A similar approach to characterizing drag in 
systems of soft grains, or in which a more viscous interstitial fluid is introduced would 
provide insight that might prove invaluable in application. 
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